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Suggested by Scullard's recent star-triangle relation for correlated bond systems, we propose a 
general "cell/dual-cell" transformation, which allows in principle an infinite variety of lattices with 
exact percolation thresholds to be generated. We directly verify Scullard's new site percolation 
thresholds, and derive the bond thresholds for his "martini" lattice (p c = l/y/2) and the "A" lattice 
(pc = 0.625457 . . ., solution to p° — 4p 4 + 3p 3 + 2p 2 — 1 = 0). We also present a precise Monte-Carlo 
test of the site threshold for the "A" lattice. 

PACS numbers: 



I. INTRODUCTION 

Percolation refers to the process of the formation of long-range connectivity in random systems 0. It has wide- 
ranging applications to problems in physics and engineering, including conductivity and magnetism in random systems, 
fluid flow in porous media, epidemics and clusters in complex networks, and gelation in polymer systems. To study 
this phenomenon, one typically models the network by a regular lattice made random by independently making sites 
or bonds occupied with a probability p. At a critical threshold p c , for a given lattice and percolation type (site, bond), 
percolation takes place. Finding that threshold exactly or numerically to high precision is essential to studying the 
percolation problem on a particular lattice. 

There are in fact relatively few lattices where p c is known exactly, all in two dimensions, where a powerful duality- 
crossing argument can be used. These lattices include bond percolation on the square lattice and site percolation on 
the triangular lattice (both with p c = 1/2), and a family of related lattices (bond-triangular, bond-honeycomb, site- 
kagome, site-(3,12 2 )) whose non-trivial thresholds can be found though the star-triangle transformation introduced 
by Sykes and Essam [2j. This transformation was successfully applied to one other lattice (the "bowtie" lattice, bond 
percolation) by Wierman 0. 

Very recently Scullard has shown that the star-triangle relation can be generalized to situations where the bonds in 
the triangular unit have correlations among them |4| . He argued that a specific correlated system could be represented 
by uncorrelated site percolation on certain lattices, leading to exact percolation thresholds for three new lattices: his 
so-called "martini" lattice, and two related lattices, which we call "A" and "B" . Scullard's proposed thresholds are 
p c (site) = 0.764826... (the solution to p 4 — 3p 3 + 1 = 0), 1/V2, and (V5 — l)/2, respectively. These rather startling 
results expand in a significant way the number and types of lattices where exact thresholds can be found. 

In this paper, we show that Scullard's arguments can be viewed in a more general way, in terms of the net 
correlation between the three vertices of a basic triangular cell. This correlation can be created by an arbitrary 
network of independent or correlated bonds, including groups of correlated bonds that effectively represent sites in 
the system. In this way, we derive the bond thresholds for the lattices Scullard considered (he already realized that 
the self-dual "B" lattice has a threshold of 1/2), give an alternative derivation of his site thresholds, and provide a 
prescription to generate additional lattices where the threshold can be found exactly. Finally, we report on numerical 
simulations to test Scullard's prediction of p c (site) for the "A" lattice. 

II. GENERAL CRITERION FOR PERCOLATION 



We consider a lattice that can be decomposed into a regular triangular array of identical triangular cells as shown 
in Fig. la, where the shaded triangles represent any network of bonds, perhaps with correlations, that connect the 
three endpoints. Randomly occupied sites can also be included within the cell by incorporating triangles of correlated 
bonds. That that array is triangular as in Fig. la is not absolutely necessary (it need only be self-dual), but we will 
consider only the triangular array in this paper. The basic triangular cell has vertices A, B, and C, and we define 
P(A, B, C) = the probability that A, B and C are all connected, P(A, B, C) — the probablity that only A and B are 
connected, P(A, B, C) = the probability that none of the three are connected, etc. 
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FIG. 1: (a) Decomposition of lattice into cells (shaded), (b) the dual-cell transformation. 



We find that the general criterion for criticality is simply 



P(A,B,C) = P{A,B,C) 



(1) 



for an arbitrary triangular cell. This formula is analogous to one that is obtained for the triangular and honeycomb 
lattices via the star-triangle transformation, but generalized to apply to triangular cells with any number of bonds 
(including correlations), arranged in a self-dual pattern such as that in Fig. la. 

We derive this result using a generalized "cell/dual-cell" transformation. The development parallels that for the 
regular star-triangle transformation, but kept in a more general context. The dual system to Fig. la is also a regular 
triangular array of non-overlapping triangles (rotated by 180°). Call the vertices A', B' , and C on the dual triangle, 
as shown in Fig. lb, and let P'(A', B' , C) etc. represent the probabilities of connecting vertices on the dual triangular 
cells. 

For the two systems to be identical, rotate the dual system so that A = A' , B = B' and C = C . Evidently, the 
two systems will have the same connectivity between vertices if we have simultaneously 



Note that for the two systems to have the same connectivity, the underlying bond occupancy of the dual cell will be 
different than that of the original cell (in fact, l—pi for a system of independent bonds), but we don't need to specify 
that in detail here. 

First we note that the two-point connectivity relations |J3J are automatically satisfied as a result of duality. We see 
from Fig. lb that if A and B connect on the original triangular cell, then A' and B' must also connect on the dual 
cell. Likewise, duality also implies that P{A,B,C) = P'(A,B,C), and P(A,B,C) = P'(A,B,C). Then, from Eqs. 
(0) and (J2J, the general condition (1) follows. This indeed represents the critical point because the two systems have 
identical connectivities, and one is the dual of the other. 

Eq. JU is equivalent to the relation P(v, h) + P(v, I, h) + P(h, l,v) — P{v,J, h) = given by Scullard for a system 
of three correlated bonds, where v, h, and I are the vertical, horizontal, and diagonal edges of the basic triangle 
drawn as a right triangle. From the identity P(v, h) = P(v, h, I) + P(v, h,l), and associating P(A, B, C) = P(v, h, I) + 
P{v, h,l) + P(v, I, h) + P(h, l,v) and P(A, B, B) = P(v,l, h), it follows that is satisfied. Note that in the present 
work we consider systems with correlations more general than those produced by three correlated bonds. 

We illustrate explicitly our argument using the usual star-triangle transformation shown in Fig. 2. Here we will 
assume that the bond occupation probabilities are pi , P2 and pz on the star (the original cell) , and P%, P2, and P3 on 
the triangle (the dual cell). Then we have explicitly for Eqs. 



P(A,B,C) 
P(A,B,C) 
P(A,B,C) 



P'(A,B,C) 
P'(A,B,C) (x3) 
P'(A,B,C) 



(2) 
(3) 
(4) 



P1P2P3 = P1P2P3 + P1P2Q3 + P1P3Q2 + P2P3Q1 

PlP2<?3 = Q1Q2P3 (><3) 

<?l92g3 + Q1P2P3 + Q2P1P3 + 13P\P2 = Q1Q2Q3 



(5) 

(6) 
(7) 
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FIG. 2: Star-triangle transformation. Left to right: star, triangle, duality orientation. 




FIG. 3: Cell, dual-cell, duality construction, and representation for Scullard's "martini" lattice (bond percolation). 



where qi — \ — pi and Qi = 1 — Pi. The two-point relation JBJ is identically satisfied if Pi = 1 — Pi- Then, there 
follows from either JSJ or Q the condition p\P2Pz — P1P2 — P1P3 — P2P3 + 1 = 0, which for p\ = P2 = P3 = P gives 
p 3 ~ 3p 2 + 1 = or p c = 1 — 2 sin7r/18 = 0.652704 ... for bond percolation on the honeycomb lattice 0]. 

Above we repeated the complete argument to find p c , working out probabilities on the dual as well as the original 
system. But using (QJ, we could just as well equate P(A, B, C) — P1P2P3 to P(A, B, C) = (719293 + Q1P2P3 + I2P1P3 + 
q_3P\P2 to find the same result without resorting explicitly to the dual lattice. 



III. THE MARTINI LATTICE (BOND PERCOLATION) 



Now we consider some new systems. First we consider the unit cell shown in Fig. 3, which is a triangle within a 
star. As shown in that figure, this cell forms precisely Scullard's martini lattice, which is a honeycomb lattice with 
triangles inserted in every other site. Here we are considering bond percolation on that lattice. We will assume that 
the outer three bonds are each occupied with an equal probability p\ , and the inner three bonds are each occupied 
with probability P2 ■ Then a straightforward calculation shows that 

P(A,B,C) = pl(p 3 2 +3plq 2 ) (8) 

P(A 7 B,C) = p\qi{p2 + q2pl)+p\{p2ql) (9) 

P(A,B,C) = q 3 + 3q 2 1 p 1 +3q 1 p 2 1 (q 3 2 +2p2q 2 2 )+p 3 1 (ql) (10) 

Equating © and (JTDJ, we find 

2pIp 3 2 - 3p\p\ - 3p\p\ + 3p\p\ + 3p\ V 2 -1 = (11) 

When pi = 1, this corresponds to a triangular lattice, and when p 2 = 1, a honeycomb lattice. When pi = P2 = P, we 
have (2p 2 — l)(p 4 — 3p 3 + 2p 2 + 1) = which yields p c — l/y/2, a new result. Interestingly, this is the same value as 
the site percolation threshold Scullard found for the "A" lattice discussed below (which is not the covering lattice of 
the martini lattice). 
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FIG. 4: Cells and construction for bond percolation on the "A" lattice. 



It is instructive to also work this out explicitly using the dual cell, which is a star within a triangle as shown in 
Fig. 3. Here Pi is the occupancy of the outer bonds and P 2 is that of the inner bonds. We find 

P'(A,B,C) = P 3 + 3P 2 Q 1 +3P 1 Q 2 (P 3 + 2P 2 Q 2 ) + Ql{Pi) (12) 
P'(A-B,C) = PiQl(Q2 + P2Ql) + Q\(PiQ2) (13) 
P'(A,B,C) = Ql(Q 3 2 + 3Q 2 P 2 ) (14) 

Then P(A,B,C) = P'(A,B,C) if Q t = p t , and equating either P(A, B, C) = P'(A,B,C) or P{A,B,C) = P'(A,B,C) 
leads to the threshold given by (|TT1) . 

Now, we can see that the transformation shown in Fig. 3 of the star-within-a-triangle to a triangle-within-a-star is 
analogous to the star-triangle transformation, but intrinsically different in that that transformation of Fig. 3 cannot be 
accomplished by applying the usual star-triangle transformation. This is an example of a more general cell/dual-cell 
transformation that can be applied to percolation problems, where the cell can be any graph connecting three vertices, 
and perhaps containing correlated bonds. 



IV. THE "A" AND "B" LATTICES (BOND PERCOLATION) 

In the same way, we can find p c (bond) for the other two lattices considered by Scullard. In Fig. 4 we show what we 
call the "A" lattice, since the basic cell has the shape of an A. It is equivalent to the cell of the martini lattice with 
the upper bond removed or made occupied with probability 1. For simplicity, we assume all bonds are occupied with 
equal probablity p. We find 

P(A,B,C) = p 2 (p 3 + 3 P 2 q) (15) 
P(A, B, C) = p 2 (q 3 ) + 2pq(q 3 + 2pq 2 ) + q 2 (16) 

Equating the two, we find p 5 — 4p 4 + 3p 3 + 2p 2 — 1 = 0, whose solution gives p c (bond) = 0.625457 . . .. This threshold 
is also new. 

Fig. 4 also shows the dual cell — a triangle with a vertical line in it — and once again, this defines a new cell/dual- 
cell transformation, and one can verify directly that P(A, B, C) — P'(A, B, C) with P — 1 — p, and rederive the above 
result for p c . 

Scullard's third lattice is generated by the cell shown in Fig. 5, which corresponds to the martini generator with two 
bonds removed. We call this the "B" lattice, because the appearance of B's when rotated by 90°, and also because it 
follows the "A" lattice. Because this lattice is self-dual, its p c (bond) equals 1/2 as noted by Scullard. This result is 
also borne out by QJ, since here 



P{A,B,C) = P (p 3 + 3p 2 q) 
P(A,B,C) = q(q 3 + 3q 2 p) 



(17) 
(18) 
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FIG. 5: Construction for bond percolation on the "B" lattice. 




FIG. 6: Construction for site percolation on the martini lattice, which is the covering lattice shown in the lower part of the 
figure (color online). 



Equating the two, one indeed finds (2p — l)(p 2 — p — 1) = 0, implying p c = 1/2. 



SITE PERCOLATION 



Now, we turn to some site percolation models. In general, some (but not all) site percolation lattices can be created 
by constructing the covering lattice of a bond problem. Here we extend this procedure by also allowing for some of the 
bonds to be correlated. Thus in Fig. 6 we show a unit cell similar to that of Fig. 3 except that now the three central 
bonds are all simultaneously occupied with probability p, and all vacant otherwise. This is indicated by coloring the 
triangle enclosed by those bonds. Note that this triangle can also be thought of as simply an independently occupied 
site, showing that this is essentially a site-bond problem, in which context this problem was first solved by Kondor 
Related correlated site-bond problems were also considered by Hu Q and Wu 0- In a somewhat different context, 
Kunz and Wu also used the idea of creating site percolation by a group of correlated bonds in bond percolation system 

El. 
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FIG. 7: Construction for site percolation on the "A" lattice (color online). 



Let pi be the occupancy of the three outer bonds, and P2 be the simultaneous occupancy of all three central bonds. 
The covering lattice (which corresponds to the equivalent pure site system) is exactly the martini lattice, as shown in 
Fig. 6. We find 

P(A 7 B,C) = p\p 2 (19) 
P(A,B,C) = qa+P2(ql+3pi<£) (20) 

and setting these equal yields (p\ — 3p\)p2 + 1 = as given by Kondor Setting p\ = p 2 = p yields p 4 — 3p 3 + 1 = 
or p c — 0.764826 ... for site percolation on the martini lattice, as found by Scullard. Here, the dual lattice is a triangle 
with a correlated star in the center, and one can verify directly that the two-point functions on the lattice and dual 
lattice are equal when Pj = 1 — pi. 

For site percolation on the "A" lattice the cell/dual cell combination is shown in Fig. 7. Here we have 

P(A,B,C) = p\ V 2 (21) 
P(A,B,C) = q 2 +P2<ll (22) 

which yields 2p\p2 — 1 = 0. When pi = p2 = p this yields p c — l/y/2 as found by Scullard. For site percolation on 
the "B" lattice, as shown in Fig. 8, we have P(A, B, C) — p\P2 and P(A, B, C) = q2 yielding (pi + l)p2 — 1 = 0. For 
Pi = P2 = P, this yields p c = (\fE — l)/2 = 0.618034 .... as found by Scullard. 

Likewise, one can create other cells that generate additional lattices where p c can be found exactly. By making the 
cell simply a shaded triangle (all three bonds occupied together), we create the site-triangular problem. A cell of just 
two bonds creates the bond-square problem. A cell with three correlated triangles touching creates the site-kagome 
problem. In this way, all known thresholds can be easily calculated, and an infinite number of new ones can be 
created, with additional lattices tending to be more and more intricate. Unfortunately, the method does not appear 
to work for some of the more notorious unsolved systems: site percolation on the square and honeycomb lattices, and 
bond percolation on the kagome lattice. 

An example of an intrinsically correlated system is one in which the cell is a simple triangle of bonds but where 
we require that in each cell at least one bond is occupied, so that P(A, B, C) = 0. Eq. implies that at criticality 
P(A, B, C) = also, so that the only non-zero correlation is P(A, B, C) and permutations, meaning that the critical 
point corresponds to all cells having exactly one occupied bond (in either a random or biased location), yielding 
p c = 1/3. If only two of the three bonds can be occupied, then the threshold corresponds to exactly one of those two 
being occupied (p c — 1/2), and we find a Scheidegger's river network model ||, many of whose properties have been 
studied and solved (e.g., fUEl)- These are examples of fixed number or canonical percolation models. 



VI. NUMERICAL TEST FOR SITE PERCOLATION ON THE "A" LATTICE" 

Finally, we report the results of a numerical test of p c (site) for the "A" lattice, whose threshold was predicted to be 
1/V2 = 0.70710678 ... 0]. We used the hull-gradient method [3[l3, with the lattice represented on a square lattice 
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FIG. 8: Construction for site percolation on the "B" lattice (color online). 
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FIG. 9: Hull-gradient simulation results for site percolation on the "A" lattice, confirming Scullard's prediction p c — l/y2 w 
0.707107. Error bars are smaller than the symbol size. 



as 2 x 2 squares with a diagonal in the left-lower corner. We considered systems of width L = 128, 256, . . . 16384 
with the gradient = 1/L. The estimates of p c (equal to the ratio of occupied to total sites in the hull) are plotted 
in Fig. El vs. 1/L, with 10 12 - 10 13 sites generated for each L. For L = 16384, we find p c (est) = 0.7071090, and the 
data extrapolate to the expected value as L — > oo with a statistical error of about ±0.000001. Note that there are 
rather unusual finite-size effects for this system (perhaps related to the orientation/representation of the lattice) as 
the estimates first increase and then decrease as L increases. 
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We did not check any of the other predictions numerically: because the cell/dual-cell argument is so compelling, 
there seems to be little doubt that these results are correct. 



VII. CONCLUSIONS AND DISCUSSION 

We have shown that the well-known star-triangle transformation usually applied to honeycomb/triangular lattices is 
a special case of a generalized transformation where triangular cells in a system are individually replaced by their duals. 
When the triangles themselves form a self-dual system as in Fig. la, then the generalized criterion for percolation 
(|T|) follows. Interestingly, to apply this criterion one does not need to explicitly apply the star-triangle or duality 
transformation, as the criterion involves connectivity probabilities on just the original system. 

Note that Fig. la is but one arrangement of the triangular cells that is self dual. Other arrangements can also be 
devised, and these will be pursued in future work. 

The basic triangular cells can possess correlations in their connectivities, and we have shown that correlations can 
be the result of having an uncorrelated bond system that is more complicated than a simple triangle or star. Site 
percolation can also be included by adding triangles of correlated bonds. In this way, we have provided a direct 
proof of Scullard's new results for the site percolation threshold for the martini lattice, and its two contractions, the 
"A" and "B" lattices. We derived the bond thresholds of these lattices also, which represents two additional new 
thresholds. Note that these three lattices are the only known cases, besides the triangular lattice, where both site and 
bond percolation thresholds can be found exactly. Other soluble lattices can easily be generated by putting additional 
bonds on the basic cells, although the lattices thus formed take on more and more the look of "decorated" rather 
than regular lattices. 

It turns out that the criterion that we derived, Eq. is closely related to work done many years ago on the 
Potts model in correlated systems 01- I n fact, Q is a special case of a criterion for being at the critical point of the 
correlated Potts model given by Wu and co-workers [ill [la. fl4l | . However, evidently the percolation limit of that 
work was never explored. 

Finally, after this paper was submitted for publication, we came across a recent preprint of work by Chayes and 
Lei in which the generalized criterion was also derived. The focus of their work is quite different, and these 
authors do not use that criterion to find the thresholds of new lattices as we have done here. 
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